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Abstract 

Our original results refer to multivariate recurrences: discrete mul¬ 
titime diagonal recurrence, bivariate recurrence, trivariate recurrence, 
solutions tailored to particular situations, second order multivariate 
recurrences, characteristic equation, and multivariate diagonal recur¬ 
rences of superior order. We find the solutions, we clarify the struc¬ 
tural background and provides short, conceptual proofs. The original 
results include a new point of view on discrete minimal submanifolds. 
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1 Discrete multitime recurrences 

The theory of multi-variate recurrences is a current effervescent topic in 
mathematics today. These recurrences are based on multiple sequences and 
come from areas like analysis of algorithms, computational biology, informa¬ 
tion theory, queueing theory, filters theory, statistical physics etc. 

We consider the lattice of points with integer positive coordinates in M n . 
A multi-variate recurrence is a set of rules which transfer a point into another, 
together with initial conditions, capable to cover the hole lattice. 
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A linear multivariate recurrence with polynomial coefficients corresponds 
to a linear PDE. In addition, extending the division to the context of differ¬ 
ential operators, the case of recurrences with polynomial coefficients can be 
treated in an analogous way. 

Bousquet-Melou and Petkovsek [2] analyse the multivariate linear recur¬ 
rences with constant coefficients (see also 0). 

Analyzing linear image processing (representations of filters), Roesser [IB] 
used a class of linear dynamical systems in two discrete-time variables, 

x{i + l,j) = A 1 x(i,j) + A 2 y(i,j), y(i,j + 1) = A 3 x(i,j) + A 4 y(i,j), 

which can be extended as block matrix system 


(Xi(t l + 1 ,t 2 , 


/ dlii(f) A 12 (t) ... Auv(t)\ 


( x i(f) \ 

x 2 (t 1 ,t 2 + 1, ...,t m ) 

= 

A 2 i(t) A 22 (t ) ... A 2 jv(f) 


x 2 (t) 

yzjv(iV 2 , + 1)/ 


\Ajsri{t) Am(t ) ... Ajvjv(^) ) 


\ x N (t) J 


where x(t) = (xi(t), ..., xjv(t)). 

Fornasini and Marchesini 0-i introduced another class of linear dynam¬ 
ical systems in two discrete-time variables, 

x(i + 1 , j + 1 ) = Ax(i + 1 ,j) + Bx(i,j + 1 ) + Cx(i,j ), 
which can be extended as 

m 

x(t + 1) = B 0 x(t) + ^ B a x(t + la). 

a=l 

This model is used for image processing, representation of discretized partial 
differential equations, models of different physical phenomena, single-carriage 
way traffic flow, and river pollution. Many other authors [I], [5], [11] , [12]-[15] 
develop the theory of filters. Both of the foregoing extensions are connected 
to mD filters theory. 

Some interesting related works are those of Prepclita 0 , ra. where a 
multiple hybrid Laplace and z type transformation was introduced to solve 
multiple differential-difference and multiple integral equations and to obtain 
the frequency-domain representations of multidimensional hybrid control sys¬ 
tems. 

The visual inspection of (i) bivariate or (ii) trivariate recurrence plots 
reveals some typical geometrical structures: (i) single dots, diagonal lines as 
well as vertical and horizontal lines (the combination of vertical and horizon¬ 
tal lines plainly forms rectangular clusters of recurrence points); (ii) all the 
above and planes. 
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The papers g], P] presents algorithms to compute stable discrete mini¬ 
mal surfaces. We add a new point of view (coming from [27], (25]) in studying 
this subject. 

In the class of multivariate sequences x{t) = xft 1 ,...,t m ) we distinguish: 
Separable multivariate sequences (those which can be written like a 
product): x(t) = x 1 (t 1 ) • • ■ x m (t m ). 

Multi-periodic sequences: 

x(t\ ..., t m ) = x(t 1 + hT\ .., t m + k m T m ), V(Aq,..., fc m ) G Z m , 

equivalently to the fact that every point (0,..., 0, tfi 0,..., 0) is a period. 

Diagonal-periodic sequences: 

x(t\ ..., t m ) = x(t 1 + T\ ..., t m + T m ), 
where T = (T 1 , ...,T m ) is a vector-period. 

2 Discrete multitime diagonal recurrence 

Any element t — (t 1 ,... ,t m ) G is called discrete multitime. A function 
of the type x : N m —» M n is called multivariate sequence. 

Let F: xlM M n and 1 = (1,..., 1) G N m . We shall study multi¬ 

variate diagonal recurrences or (discrete multitime diagonal finite difference 
equations ) of first order 


x(t+ i) = F(t,x(t)), (1) 

where the multivariate vector sequence x(t) = (x 1 (t),..., x n (t)) is solution 
of the system (jT|) . 

This model of diagonal recurrence can be justified by the fact that to a 
first order PDE system 


Dix{t) = f(t,x(t)), t G M m 
we can associate a discretized equation of the form 

x(t + 1) — x(t) = f{t , x(t)), t G Z' m . 

The initial (Cauchy) conditions on a curve, for the PDE system, are trans¬ 
lated into initial conditions for the diagonal recurrence. 
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Remark 1. Let us consider a two-variate recurrence with the unknown se¬ 
quence x(m,n). A very frequent and interesting case happens when the dif¬ 
ference m — n, between the arguments of the unknown sequence, is constant 
among all its occurrences in the multivariate recurrence relation. For in¬ 
stance, this happens for any diagonal recurrence x(m, n) = f(x(m — 1 , n — 1 )) , 
where the difference between the first and second argument of x(m, ri) is al¬ 
ways m — n. Such a recurrence can be rewritten as a univariate recurrence 
y(t ) = f(y(t — 1)) ; where y(t — k) = x(m — k,n — fc), V 7 G N. Another inter¬ 
esting case, similar to the one above, is when the sum m + n of the arguments 
of the unknown sequence is constant. For instance, multivariate recurrences 
of the form x(m,n ) = f(x(m + 1 ,n — 1 )) can be rewritten as univariate 
recurrences y(t ) = f(y(t — 1)), where y(t — k) = x(m + k,n — k),Wk G N. 

2.1 Linear discrete single-time recurrence 

Let us recall a well-known result regarding a single-time linear recurrence 
equation, i.e., m = 1 , in an original version that can be extended to the 
multi-temporal case. 

Proposition 1. Let A: N —>■ At n (K), b: N —* M” and x 0 G Then the 
unique sequence x: N —> M n which verifies the first order linear single-time 
recurrence equation 


x(t + 1) = A{t)x{t) + b(t), Vt G N, 

x(0) = Xq 


( 2 ) 


is 


x(t) = Aft — 1 )A(t — 2 ) • ... • d(l)yl(0)xo+ 

t—2 

+b(t - 1 ) + ^2A(t- l)d(t -2 A(k + 1 )b{k), Wt > 2, (3) 

k=0 

x(l) = d( 0 )x 0 + b( 0 ). 

Proof. Mathematical induction after t > 1. For t — 1, the result is verihed 
automatically, using the hrst relation in (| 2 ] 1 , for t = 0 , as well as x(0) = xq. 
For t — 2, x(2) = d(l)x(l) + 6(1) = d(l)d(0)o:o + d(l)6(0) + 6(1). 

It remains to solve the inductive step: prove that, if the statement holds 
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for some natural number t > 2 , then the statement holds for t + 1 ; indeed, 

x(t + 1 ) = A(t)x(t) + b(t ) 

=A(t)A(t - 1 )A(t - 2) • ... • y4(l)A(0)x 0 + 

t—2 

+A(t)b(t — 1 ) + ^2 A(t)A(t — l)A(t — 2) • ... • A(k + 1 )b(k) + b(t) 

k =0 

t -1 

—A(t)A(t — 1 ) •... • A(l)v4(0)xo + 6 (f) + — 1 ) •... • d(A’ + l)fc(fc). 

k =o 


□ 

Corollary. Let A G .M n (R) 6e a constant matrix, and xq G M n . Then the 
unique sequence x: N —> R n which verifies 

x(t + 1) = Ax{t) + b(t), Vt G N, 
x(0) = so, 

is 

t- i 

x(t) — ^t-i-k 5 ^^ Vt > 1. (5) 

k=0 


2.2 Linear discrete multitime diagonal recurrence 

The diagonal discretization of PDEs incorporate points which he on diagonals 
of the grid. This is a good enough reason for introducing and analyzing the 
diagonal recurrences. 

Let m > 2 , A: —y At n (R), b: N rn —> R n . In this Subsection we refer 

to a linear discrete multitime diagonal recurrence equation of first order 


x{t + 1) = A(t)x(t) + b(t), t G N m , ( 6 ) 


with x : N m —* R n = At n ,i(R). 

For convenience, we denote //(t) = minjt 1 , t 2 ,..., t m }. 


Lemma 1. Let D 


t = (t\t 2 ,...,t m ) G N m 


n{t) = t m 


m > 2 , and 


A: D —>• At n (R), b: D —>■ R n , /: N TO_1 —* R n . Then the unique multivariate 
sequence x : D —> R n which verifies the first order linear discrete multitime 
diagonal recurrence equation 


x(t + 1 ) = A(t)x(t) + b(t), 

x(t 1 ,...,r- 1 ,o) = /(t 1 ,...,r- 1 ), 
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Vt G D, 

V^ 1 ,...,^- 1 ) G N m_1 


(7) 



is 

x(t) = Aft - l)A(t - 2 • 1 ) • ... -A(t - t m • 1 )f(t 1 - t m ,..., f 71 ' 1 - t m ) 

t m 

+b(t — 1 ) + H(f — l)H(f — 2 - 1 ) • ... • A(t — (k — 1 ) • 1 ) 6 (£ — k ■ 1), ( 8 ) 

k=2 

Vt G D, with t m > 2, 

..., r- 1 , i) = - i, ..., r - 1 - i, o)/(t x - i, , r - 1 - i) 

+ 6 (t 1 - l ,...,*"* -1 - 1,0), V(f\..., £ m_1 ) G (N*) m_1 . U 

Proof. Let us remark that for any t — ft 1 , t 2 ,..., t m ) G D and any s G N, it 
follows t — ( t m — s) ■ 1 = (t 1 — t m + s,..., f m_1 — t m + s, s) G D. Hence, for 
a fixed t E D, we can define the sequences 


H:N^ 2 W n (R), 

x: N -)• R n , 


l(s) = H(£ - (t m - s) • 1 ) = ^(f 1 -t m + s,. 

..,r-i -t m + s, s ), 

( 10 ) 

b(s) = b{t - {it m - s) • 1 ) = bit 1 -t m + s,.. 

-t m + S,s), 

( 11 ) 

x{s) = x(t — (■ t m — s) ■ 1 ) = x{t 1 — t m + s,. . 
The sequence x verihes 

,.,r-i -r + s,s). 

( 12 ) 

x(s + 1 ) = H(s)x(s) + b(s), 

Vs G N, 

(13) 

x( 0 ) = f{t 1 - t m , ...,t m ~ 



Now we can apply the Proposition [lj For any s > 2, we have 
x(s) = A(s — l)H(s — 2 ) • ... • A(0)x(0) + 

s -2 

+b(s — 1) + H(s — l)H(s — 2 ) •... • A(k + 1 )b(k) 
k =o 

i.e., 

x(t-(t m -s)- 1) = A(t-(t m -s+l)-l)A(t-(t m -s+2)-l)-.. .■Aft-t m - l)x(0)+ 

+b(t -{t m -s + 1 ) • 1 ) 

s—2 

k=0 
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For t m > 2, we set s = t m in the foregoing relation. Then 

x(t) = A(t — 1 )A(t — 2 • 1 ) • ... • A(t — t m ■ l)x(O) + b(t — 1 ) 

t m -2 

+ 1W - 2 • 1 ) •... • A{t - {t m - k - 1 ) • 1 )b(t - ( t m - k ) • 1 ) 

k=0 

= A(t - 1 )A(t - 2 • 1 ) •... • A{t - t m • 1 )f(t 1 - t m ,..., t™- 1 - t m ) + b{t - 1 ) 

t m 

+ J2 A (t~ 1 )A(t - 2 • 1 ) •... • A(t - (t m - k + 1 ) • 1 )b{t - ( t m - k + 2 ) ■ 1 ) 

k=2 

= A{t - l)A{t - 2 • 1 ) •... • A(t - t m • 1 )f(t 1 - f m ,..., t m ~ l - f m ) 

t m 

+b(t - 1 ) + ^ A[t - 1 )A(t - 2 • 1 ) •... • A(t - {k - 1 ) • 1 )b(t - k • 1 ) 

k =2 

For = 1, we obtain the formula (l9|) : 

x(t\..., f m_1 , 1 ) = xiit 1 - 1 ,..., t m ~ l - 1 , 0 ) + 1 ) = 

= A{t x -1,..., t m_1 -1, o)x(t 1 -1,..., t m_1 -i,o)+ bit 1 -1,..., t™- 1 -1, o) 
= Ait 1 - 1 ,..., t m_1 - i, o)/^ 1 - i, ..., r -1 - 1 ) + bit 1 - 1 ,..., r- 1 -i,o). 

□ 


Remark 2. In the conditions of Lemma [lj if moreover the function Tl(-) is 
constant, i.e., A{t) = A, Wt, then the formulas (jSJ) and (jH)) become 

t m 

x{t) = Al" 1 fit 1 - t m ,..., t ™- 1 - t m ) + ^ - k • 1 ), (14) 

k= 1 

Vt G D, with t m > 1. 

Notation: for every /3 G {1,2,..., m}, we denote 

it 1 ,...,t m ) ■= it\...,tP-\tP +1 ...,t m ) G FT -1 . 


Theorem 1. Let m > 2, A: FT" 1 —» A4 n (M), b\ N m —>• M n . VFe consider the 
(m — 1 )-sequences fi, / 2 ,..., f m : srtc/i that 




Vt 


i 

5 


y.a — 1 jlQ+1 
* * 1 1 I 1 1 


t a = 

...,t 0 ~\t 0+1 ,...,t m G N, 


(15) 
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for any a, (3 G {1,2, Then the unique m-sequence x: —>■ M n 

which verifies 


x(t + 1) = A(t)x(t ) + b(t), Vt G N m , 

x(t) = fpty 1 ,...,t m ), V(f\ 
tP=o 

VP G {1,2,..., m}, 


e N 


m —1 


is defined either by the formula 

x(t) = A(t - 1 )A(t - 2 • 1 ) • ... • A(t - t p • 1 )- 

• fyif 1 -t/ 3 ,. t m ~ l - t p ) + b(t - 1 ) 

tP 

+ ^2 A{t — 1 )A(t - 2 • 1 ) •... • A(t - (k - 1 ) • 1 )b(t -k- 1 ), 

k =2 

if fi(t) =tP> 2 , 


or by the formula 

x(t) = A(t' - 1 1 - l,0,i' 3+1 - - 1). 

■M? 

+ b(t‘ - 1,..., t 3-1 -1,0, t 3 * 1 - 1,1), 
if rft) = t 0 = 1 . 


(16) 


(17) 


(18) 


Proof. Let us remark that the multivariate sequence x is well defined, i.e., 
if t a = t 13 = /i(t), then the expressions which define x of the formulas (TT7T) 
(or (USD), corresponding to a and coincide due to the equalities (THU) . 

If t 13 = n(t), the conclusion follows applying directly the Lemma [TJ having 
instead of t m . 

Conversely, one observes immediately that if the multivariate sequence x 
is defined by the formulas (fTT|) . (fT 8 l) . together with the second formula of (1T6|) 
(for tfi = 0 ), then the multivariate sequence x verifies also the first relation 
in (USD , i.e., 

x(t + 1) = A(t)x(t) + b(t), Wt G N m . 


□ 


Remark 3. Let m > 2, A: N m —>■ _M„(R), b: N m —> M n and the (m — 
1 )-sequences /i, / 2 ,..., f m ■ M m_1 —> M n . If the m-sequence x: N m —>• M n 
verifies the relations 


x(t) V(( 1 ,...,« 3 ,...,r) er-‘, 

tfi= 0 

V/3 6 {1,2,...,m}, 


( 19 ) 




then, 


up,...,n 

tP =o 


t a =o 




( 20 ) 


This follows immediately since, x(t) 


tP=0,t a =0 


x(t) 


t a =0,tP=0 


Corollary. In the conditions in Theorem [H if moreover, the function A(-) 
is constant, i.e., Aft ) = A, \/t, then the formulas (fT7|) and (jT8l) become 


x(t) = A tf> fpft 1 

if fi(t) = tP> 1. 


tP 

t?) + ^A k ~ l bft - k ■ 1 ), 

k =1 


( 21 ) 


Identifying the initial conditions with the constant diagonal recurrence, 
we obtain new information about the discrete diagonal flow. 

Proposition 2. Let m> 2 and A: N m —» Ad n (M). Denote 


S 


y: N m -> 


y(t + 1) = j/(t),Vt G N n 
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y = -)• 

and introduce the function 


x(t + 1) = A(t)x(t),Vt G N" 

ip: S^V, 





■y(t- y(t) • 1), 


if > 1, 
*/ M*) = o. 


( 22 ) 


a) T/ie sets S and V are real vector spaces, and if is an isomorphism of 
vector spaces. 

b ) The vector space V has infinite dimension. 


Proof. a) First, let us observe that the application f) is well defined, i.e., 
the m-sequences definite by the formula (|22|) verify the recurrence x(t +1 ) = 
A(t)x(t). This follows immediately from the Theorem [0 

The respective sequence is the unique sequence x : N m —> M n which veri¬ 
fies the problem 


x(t + 1) = A(t)x(t), Vt G M m , 


x(t) 


= y(t 1 ,...,tf 3 \0 ,t m ,. 

tP =o 

V/3 G {1,2,..., m}. 


,r), v(tV.., * m ) gn 


m —1 


? 


(23) 
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Obviously, V and S are vector spaces over R. 

One observes that 0 is a morphism of vector spaces. Let us prove that 
the function 0 is injective: let y E S such that 0(t/(-))(f) = 0, Vf E LT m . 

ft follows 0 (y(-))(f) = 0 ; but, according to the relation (l? 2 j) . we have 


tfi= o 

= y(t) and we obtain y(t) 
v ' tP= o tP= o 

m-sequence which verihes the problem 

y(t+l) = y(t), Vte N m 

y(t) = o, v(t\...,t0. 

tP=0 

V0 E {1,2,..., mj. 


tP =o 


= 0. Hence y is the unique 


, t m ) e N' 


m—1 


Applying the Theorem [H it follows y(t) = 0, Vf E N m , i.e., y is m-sequence 
zero. 

The surjectivity of 0: let x E V. We choose y: N m —* M n , the unique 
m-sequence which verihes 


y(t+ 1 ) = y(t), 


y(t) 


= x(t\. 

tP =0 

V/3 E {1,2, 


Vf E N m , 

.y-\oy + \ 

. ,,m}. 


,f), V^ 1 ,...,^,...,^) e N 


m—1 


Obviously y E S and x(t ) 


t/3=0 



Hence the m-sequence x verihes 


the relations (J23]). It follows (Theorem Q]) that 


x(t ) = A(t — 1 )A(t — 2 • 1 ) •... • A{t — L 3 • 1 ) • y(t — t 13 ■ 1 ), if /i(t) = t 13 > 1 ; 


z(f) = 2/(t), if hO) = 0, 

hence, according the relation fl22jl . we have 0(y(-))(f) = x(t), Vf E N m , i.e., 

0 (y(-)) = ®( 0 - 

6 ) Since V and S are isomorphic vector spaces, it is sufficient to show that 
S has an infinite dimension. Equivalently, we shall show that S contains an 
infinity of linearly independent elements. 

Let v E R", v 0 0. For each k eN*, we consider the sequence 

y k : LT R", y k {t) = ( t 1 - t 2 ) k v, Vf = (f 1 , f 2 ,..., t m ) E hP. 

Since y k (t + 1) = (f 1 + 1 — f 2 — l) fc u = t/fc(f), it follows y k E S. 

Let F C N*, F hnite and non-void. For each k E F, we consider the 
sequence a k E R, such that = 0, i.e., a k yk(t ) = 0,Vt G LP. 

fceF fceF 
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Setting t 2 = 0, it follows ^ a^t^v = 0, Vf 1 6 N, or | jn = 

k£F ' keF ' 

OjVt 1 G N. Since v ^ 0, we deduce that ^ Qfc(t 1 ) fc = 0, Vf 1 G N. 

k&F 

Consequently the polynomial P(X) := cikX k vanishes for any t 1 G N, 

keF 

i.e., P(X) has an inhnity of roots; whence, it follows that P(X) is zero 
polynomial, i.e., a*, = 0, Wk G F. 

We proved that [yk \ k G N*} is a subset of S consisting in linear inde¬ 
pendent elements. Hence the dimension of S is infinite. □ 

Remark 4. For determining the solutions of the recurrence 

x{t + 1) = A(t)x(t), t G LT m , 


it is sufficient to know the solutions of the recurrence 

y(t + l) = y(t), tGr. 

Suppose that the functions /i, / 2 ,. .., f m '■ —> HP satisfy the rela¬ 

tions (USD- We want to determine the multivariate sequence x\ N m —>■ M n 
such that 

x(t + 1) = A(t)x(t), Vf G N m , 

= /»(«', n, v(t*, netr-\ 

V/3 6 {1,2,...,m}. 

If y: N m —* W 1 verihes 

y(t + 1) = y(t), Vie N m , 

y(t) =Mt\.n, v(i 1 ,...,iV..,r)epr- 1 , 

tP =0 

V/3 G {1,2,...,m}, 
then, from the proof of the Proposition O it follows x(-) = 


2.3 Examples 

Bivariate recurrences Let A G Al n (M) and two sequences /i, f 2 '■ N — > M n , 
with /i(0) = fi (0). Then the unique double sequence x: N 2 —> which 
solves the problem 

( x(t 1 + 1, t 2 + 1) = A ■ xft 1 ,t 2 ), V(t 1 , t 2 ) G N 2 , 

\ x(0 ,t 2 ) = /r(t 2 ), Vt 2 GN, 

I z(t\0) = f 2 (t 1 ), Vt 1 G N 
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IS 


xft 1 , t 2 ) 


A* 1 fi(t 2 — t 1 ), if F<t 2 
A* f^t 1 ~t 2 ), if t 2 <t\ 


Trivariate recurrences Let A G AT 


and three double sequences 


/i,/ 2 ,/ 3 :N 2 ^ 


i n , with 


/i(0,t 3 ) = /2(0,t 3 ), /i(f 2 ,0) = / 3 (0,t 2 ), / 2 (f\0) = / 3 (*\0), Wt\t 2 ,t 3 G N. 


Then the unique triple sequence x: N 3 —> W l which solves the problem 


' xft 1 + 1, t 2 + 1, t 3 + 1) = A ■ x(t l ,t 2 , t 3 ), V(t 1 , f 2 , t 3 ) G N 3 , 

x(0,t 2 ,t 3 ) = f 1 (t 2 ,t 3 ), V(t 2 ,t 3 ) e N 2 , 

z(t\0A 3 ) = f 2 (t 1 ,t 3 ), V(t\t 3 ) e N 2 , 

x(t\t 2 ,0) = f 3 (t\t 2 ), V^d 2 )^ 2 


is 


^(f 1 , t 2 , t 3 ) 


A* 1 hit 2 -t\t 3 -t 1 ), 

A t 2 f 2 {h-t 2 ,t 3 -t 2 \ 


if t 1 — min-ft 1 , t 2 ,t 3 } 
if t 2 — minjt 1 , t 2 ,t 3 } 
if t 3 = minjd 1 ,t 2 ,t 3 }. 


2.4 Solutions tailored to particular situations 

To find the solution for a discrete multitime diagonal recurrence, with con¬ 
stant coefficients, x(t + 1 ) = Ax(t), we can use a discrete single-time re¬ 
currence y(t + 1) = Ay(t ), together a family of initial conditions (see the 
foregoing examples). Particularly, for any constant matrix A, the solution of 
first order discrete multitime diagonal recurrence can be written as 

x(t) = A <£,t> x 0 , € = (ei,e m ), < e, 1 >= 1. 

The existence conditions of the powers of the matrix A gives the conditions: 
(i) e G Z m if A is non-degenerate, (ii) e G if A is degenerate. If we add 
initial condition, then, in both cases, one and only one component of e is 
non-zero (depending on the initial condition). 

Theorem 2. Let A G Af n (M) be a matrix which verify the equation A m = A. 
For any x 0 G M n (constant), the vector x(t) = A t + --- +tm x 0 verify the diagonal 
recurrence x(t + 1) = Ax(t). 

Proof. By computation, x(t + 1) = A m+tl+ - +tm x 0 = A m A tl+ --- +tm x 0 = 
Ax(t). □ 
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A matrix B is said to be an m-th root of an n x n matrix A if B m = A, 
where m is a positive integer greater than or equal to 2. If there is no 
such matrix for any integer m > 2 , then A is called a rootless matrix. A 
non-singular matrix and a diagonalizable matrix have m-th roots in complex 
numbers (see BHD- 

Theorem 3. If the n x n matrix A has an m-root B, then a solution of the 
diagonal recurrence x(t + 1 ) = Ax(t) is x(t) = B* +--- +tm x 0 , where xq is a 
constant vector. 

Proof. Explicitly, we have x(t + 1) = B m + tl +-+ tm Xq = B m B tl+ - +tm x l0 = 
Ax(t). □ 

2.5 Characteristic equation 

We refer to the recurrence 


x(t + 1 ) = Ax(t), 

satisfying the initial conditions (TIBI) made compatible by (1151) . 

If the matrix A is diagonalizable, then there exists a basis {vi,...,v n } 
consisting in eigenvectors, in C n = At ni i(C). Denote by X k the corresponding 
eigenvalues (distinct or not). Any solution of the previous recurrence is of 
the form 

n 

x(t) = YMt) v k, 

k =1 

where Cfe(t) G C and Ck(t + 1) = Ck(t). The functions C}ft) are determined 
by the initial conditions (II til) . 

Remark 5. We look for solutions of the form x(t) = v \ <e,t> , with < e, 1 >= 
1 and u / 0. The existence conditions of the powers of the eigenvalue X 
gives the conditions: (i) e G if A is non-degenerate, (ii) e G if A is 
degenerate. It follows Av = Xv. Consequently, X is an eigenvalue, and v is 
an eigenvector. The equation P(X) = det(A — XI) = 0 is called characteristic 
equation. If we can determine n pairs (A k,Vk), then a particular solution of 
the recurrence is 

n 

x (t ) = Y° k ^T’ t>V k- 

k =1 

When e = (0,..., 0,1, 0,..., 0), the scalar product is < e,t >= t a . We find 
a solution of the form 

n 

x (t) - y c k Af v k . 

k=\ 


13 


3 Multivariate diagonal recurrences of 
superior order 

Let k > 2 and the matrix sequences 

B k _ x : LT M n {R),f: LT -G ML 

The linear diagonal vectorial recurrences of order k have the form 


k -1 

x(t + k-l) = Y^ + J ■ !) + f{t), (24) 

3=0 

with xit 1 ,..., t^ _1 , 0 , t & +1 ,..., t m ), .. xit 1 ,..., h 3-1 , fc— 1 , h 3+1 ...., t m ) given 
for any /? G {1,2,..., m}. The unknown sequence is x : N m —>• M”. 

These recurrences easily be reduced to order one recurrences of the form 
(USD- Indeed, it is enough to consider 


y: r (R n ) k = AW(M), 6: -G (M”) fc = A4 nM (M), 



( 


\ 


( On, 1 \ 




x(f + 1 

) 


O n .i 


2/(f) = 



a;(i + 2 • 

1) 

; Kt) = 


5 






071,1 



\ + (fc - 1) • !) / 


V fit) / 




A: N 

m ^M 




/ 


In 

O n 

... O n 

0 n 

\ 




O n 

In 

... O n 

O n 

A(f) = 










O n 

O n 

... 

In 


V 

Biit) 

B 2 (t) 

• • • B k _ 2 (t) 

Bk-i(t) J 


Then the sequence y verifies 

y(t + l) = A(t)y(t)+b(t), Vt€fr\ 
being given yft 1 ,..., t^~ l , 0, C _l ..... /"'), V/? € {1, 2,... ,m}. 


(25) 


(26) 


(27) 
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4 Second order multivariate recurrences 


Apart from the first order multitime recurrences, an important role is played 
by the second order multitime recurrences. 


Lemma 2. Let A G _A/f 2 (M) and Ai, A 2 its eigenvalues, i.e., the roots of the 
polynomial 

P{\) = A 2 -Tr(A)A + det(A). (28) 

i) If Ai, A 2 G M and Xi ^ A 2 , then \/k G N 


A k 


\k \ k \ \k \ \k 

7 '1 ^2 a ^ 2^1 t 

~7-T - A -t- 7 - 


ii) If Ai = A 2; f/ien VA G N 


(29) 


A k = kX*- 1 A- (A- l)AjI 2 . 


(30) 


zzz) If A 1; A 2 G C \ M, Ai, A 2 = r(cosd ± z sin 9), with r > 0, 0 G (0, 2zr) \ {zr}, 
then VA G N 

r fc- 1 sinAd^ r k sin(A — 1 ) 6 * ^ 
sin 0 sin 9 

Proof. There exists Q( X) G M[AT] and cq, ci G 1 such that 

X k = P(X)Q(X) + Cl X + c 0 . 

Hence A k = P(A)Q(A) + ciA + c 0 / 2 . Since P(A) = 0 2 , we obtain 

= ciA + cq/ 2 . 


If Ai ^ A 2 , for X = X U X = X 2 ( P(Xj) = 0 ), then 

Ai = C1A1 + cq, A 2 = ciA 2 + cq, 


whence we find 

\ k \ k \ \ k \ \ k 

_ Ai — a 2 _ a 2 a x — Ai a 2 

Cl “ Ax - A 2 ’ Co “ Xi — A 2 ' 

Consequently, we get the results in i) and in). 

If Xi = A 2 , then Ai is a double root of the polynomial S(X) := X k — 
Ci A' — Cq. Hence S( Ax) = 0 and S'(Xi) = 0. From these relations it follows 
Ci = AA* -1 , c 0 = —(A — l)Ai and one obtains the result in ii). □ 
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We consider now diagonal recurrences of order two 
x(t + 2 • 1 ) + ax(t + 1 ) + bx(t) = 0 , 


( 32 ) 


with o,fteK constants and with initial conditions 

x(t l ,..., V 3-1 , 0, t p+1 ,.. •, t m ), x(t\ ..., V 3-1 , 1, t ni , 

given for any /3 G {1,2,..., m} and x : FT m —> R unknown sequence. Denoting 


y(t) = 


x(t) 

x(t + 1) 


A = 


0 1 
—b —a 


the multivariate sequence y verifies y(t + 1) = Ay it). From the Corollary of 
the Theorem [TJ it follows 


«'«=(* + i))=^»(‘-^- 1 ) = A*'( 


x(t) 
x(t + 

if fi{t) = t 13 . 


x(t -t^- 1 ) 

x(t — (V 3 — 1) • 1) J ’ 


The characteristic polynomial P(X) = A 2 + a\ + b of the matrix A has 
the roots Ai,A 2 - According the Lemma [21 the matrix A tf> is of the form 
A t,P ‘ = c\(t^)A + c 0 (t^)l 2 (formulas (j2Ul) . (j50|) . (IBTT) 1. Hence 


) A _ Cl (t P\ A ( x(t-tP-l) \ , p, ( x(t - tf* ■ 1) 

1) ) ~ Cl(t 3 A ^ x(t - (V 3 - 1) • 1) ) + 3 ^ x(t - - 1) • 1) 


x{t) 
x(t + 

if n(t) = . 

Consequently 

x(t) = Ci(t l3 )x(t — ( t 13 — 1 ) • 1 ) + c 0 (t l3 )x(t — ■ 1 ), if fj,(t) = t 13 . 


We have proved the following result 

Theorem 4. Let m > 2, a, b G R and Ai,A 2 the roots of the polynomial 
P( A) = A 2 + aA + b. Suppose that the (m — 1)-sequences 

h,h, M ™- 1 -G R, gi ,g 2 , N m_1 -G R, 

satisfy, for any a, (3 G {1,2,..., m}, the compatibility conditions 


fait 1 ,- 

fa 

• • J 1 5 * 

,.,n 

=* (<i - • • 


■ ,n 

t a = o’ 


g^t 1 , 

fa 

• * ? 1 j • 

..,n 

• ■ 


.,n 

, 

t a = 1 

(33) 

fai^,. 

fa 

• * 1 1 > • 

..,n 

■ ■ 


.,n 

t“=o’ 



Vi 1 , • 

■■,t a ~ 

1 «Q+1 i/3-l 

,6 , . . . , 6 

t^\. 

..,t m 

G N. 
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Then the unique m-sequence x : N m —* M which verifies 
x(t + 2 • 1 ) + ax(t + 1 ) + bx{t) = 0 , WtE N m , 

x(t) = ...,£,■■■>n, n G r-\ 

V =0 

X(t) = n,...,r), V(f \. t m ) e N m_1 , 

n=i 

V7 e {1,2,... ,m}, 

is defined by the following formulas: 
i) If Ai, A 2 el and X\ ^ A 2 , i/ien 

x(t) = ^gpit 1 -t p + 1 ,..., P,..., f "' 1 - t p + 1 ) 


Ai — A 2 

\ \t& \ \tP 

A 2 A x - AiA 2 , 1 


./sR -t' 


m— 1 _ #/3 


*/ M*)= 


Ai — A2 
ii) //Ai = A 2 , £/ien 

z(0 = t l3 \f-' g f (t 1 - t 0 + 1,..., P,..., f"- 1 — t' 3 + 1) 

-(t fl -i)Af /(,(*'- 1?,■ ,*” _1 -t s ), i/ MO = t§ ■ 


( 34 ) 


(35) 


(36) 


m) If A 1 ,A 2 eC\R, Ai, A 2 = r(cos 6 ± isin 0), with r > 0 , 9 e ( 0 , 2 - 7 r) \ { 7 r}, 
then 

/ x r^ _ 1 sint /3 0 „ 

— 9p(t n t p + 1,. . ., tf,. . ., t - t p + 1) 


sin 6 


r l sin (tP — 1)6 
sin 6 


(37) 


fpit 1 ~ t p ,...,..., t m 1 - t p ), if n(t) = t! 3 . 


5 Discrete minimal submanifolds 

Let (M, g) be a Riemannian manifold and (N. h ) be a Riemannian oriented 
submanifold (possibly with boundary). Let x = (ad), i = 1 be the 

local coordinates in M and t = (t a ), a = 1,..., m, the local coordinates in N. 
If the parametric equations of the submanifold N are x l = x l (t), then the 
induced Riemannian metric has the components 

Kp (t) = gij{x{t))x' a {t)x?p{t), 

where x l 0 (t) = |p)(£). The determinant of this metric is denoted by d — 
det(h a p(t)). If £ C N is compact subset, corresponding to t £ Q- compact, 
then its area is 

[ da — f Vd, dt 1 A ... A dt m . 

is Jn 
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The submanifold £ is called minimal if and only if it is a critical point of the 
area functional 

I(x(-)) = [ Vd dt 1 A ... A dt m , 

Jn 

for all compactly supported variations. Introducing the Lagrangian 

L — Vd — \J det (gij(x{t))x* Jl {t)x 1 p(t)) 
a minimal submanifold is solution of Euler-Lagrange PDEs system 


dL dL 

dx i Q dx^ 


0, 


i.e., vanishing mean curvature vector , 


r 


o. 


Though continuous models of minimal submanifolds are usually more conve¬ 
nient and yield results which are more transparent, the discrete models are 
also of interest being in fact discrete dynamical systems. 

The theory of integrators for multi-parameter Lagrangian dynamics shows 
that instead of discretization of Euler-Lagrange PDEs we must use a dis¬ 
crete Lagrangian, a discrete action, and then discrete Euler-Lagrange equa¬ 
tions (see DSD- Of course, the discrete Euler-Lagrange equations associated 
to multitime discrete Lagrangian can be solved successfully by the Newton 
method if it is convergent for a convenient step. 

To simplify, we consider the minimal 2-dimensional submanifols, having 
the coordinates t = (t 1 , t 2 ). The discretization of the Lagrangian L(x(t), x a (t)) 
can be performed by using the centroid rule (see [19], [27], [2B]) which con¬ 
sists in: (i) the substitution of the point (t 1 , t 2 ) with (mh 1 , nh 2 ), for the fixed 
step (h 1 ,/? 2 ); (ii) the substitution of the point xlpt 1 ^ 2 ) with the fraction 

>. _ X m n n %mn -\-1 _ 

S mn g i 

and (iii) the substitution of the partial velocities x\ = J-p-, x 2 = J-p- by the 
fractions Xm + 1 : i ~ rZm " , x ™+}- x ™ m We can write 

hi ’ h 2 

L 2 d = det h a p(m, n ) 


and since 


ha/3 Qij^jLmn) 


[xinn )/3 


h c 


hP 


Vmn)l ^m+ln ^ mn) (^ran)2 X r 


mn+1 


-x r 
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it follows 


h\h 2 2 L 2 d = det • 

The discrete Euler-Lagrange equations are 

E^© = ° 

^ '-'•A'mn 

where £ runs over three points, 

( K x mni x m +! n • x mn+ 1 ) 

(^m~l n.) ^m ni ^m-ln+l)i left shift map 


and 


Since 


d(x. 


dx L„ 

explicitly, we compute 

dL d 


(y^m n— 1 5 x m+in-i, x mn ), right shift map 
m — — 1, n — 1,..., N — 1. 


mnJa _ X* ^m-ln)« _ £l ^{ X mn—l)oi _ x* ^2 

~ dx* ~ fcd “’ dx* “ fc 


f)qr>k 
w mn 


(y%mm ^m+lm ^mn+l) 


1 <9d <9/i 


O'/? 


W <9 /?. a/3 dx k mn 

= Vd,h a/3 n )a«J/3 - g kj (« n )a + «„)/?)) * 


9 ^ 1 N 

" \^m —1 ni %mm %m— 1 n+ly 


i % 


= \fdh ap ( --^j:{U-in){.x l m _ ln ) a {x J m _ ln )p + g kj {(x 3 m _ ln ) a S} + S^x 3 ^ 


9 f \ 

" V'^mn—1 5 ^m+1 n— 1 5 %mn) 


— \fdh P ( l (£mn-l)(x mn _i) a (x 3 mn _ 1 )i3 + g k j ( (x J m m _i)a^/3 + ^ a ( x mn_l 


3 drc fc 


19 



Theorem 5. Denoting 


Amnafi 9 kj (( X mn)oc T ip^mnlp) 

A-m-lnap — 9kj (( X mn)a &f) + d Q ( X mn)p ) 

-A-m n—lafi = 9kj (( X rnn)a + $a { X mn)p) ) 

the variational integrator of discrete minimal 2-submanifolds is described by 
the recurrence equation 


X 

m,m—l,n—l 




lTl,n) J — k (£rnn)(x mn )a(Xm n )/3 + A mria p 


0 . 


6 Conclusions 

This paper presents original results regarding the multivariate recurrence 
equations. Our approach to multivariate recurrence equations is advanta¬ 
geous for practical problems. The original results have a great potential to 
solve problems in various areas such as ecosystem dynamics, financial model¬ 
ing, economics, image processing (representations of filters), and differential 
geometry etc. The two-dimensional filters are extensively used in processing 
two-dimensional sampled data (seismic data sections, digitized photographic 
data, gravitational and magnetic maps etc). 
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